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Abstract 

After presenting a survey of theoretical results concerning the structure 
of two-dimensional QCD, we present a numerical study related to the mass 
eigenstates and the decay amplitudes of higher mesonic states. We discuss in 
detail the fate of important dynamical points such as stability of the spectrum 
and the problem of screening versus confinement in this context. We point 
out differences in the large distance behaviour of the potential, which can be 
responsible for the question of stability of the spectrum, as well as whether it 
is finite. 
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1 Introduction 



Unlike the Schwinger model (p, Quantum Chromodynamics of massless 
fermions in two dimensions is not exactly solvable 0. It nevertheless serves as 
a very useful laboratory for studying problems such as the bound-state spec- 
trum and algebraic structure. These problems are important tools for general 
understanding of realistic quantum field theories and are expected to realize 
the important features of four dimensional Quantum Chromodynamics. In 
particular, exact properties may be derived, once one arrives at an equivalent 
bosonic formulation in the form of a gauged Wess-Zumino-Witten (WZW) 
action |, @. 

The first attempt to obtain the particle spectrum dates back to 1974, and 
was based on the 1/N expansion^, |6|, where N is the number of colours. In 
this limit one is led to a bound state spectrum corresponding asymptotically to 
a linearly rising Regge trajectory. The use of the principal-value prescription 
in dealing with the infrared divergencies is however highly ambiguous due to 
the non-commutative nature of principal-value integrals. Moreover, the result 
for the fermion propagator is tachyonic for a small fermion mass as compared 
with the coupling constant, hence, in particular, for zero fermion mass. This 
has made 't Hooft's solution a controversial issue@- 

In the large N approximation the gluons remain massless, since fermion 
loops do not contribute to the Feynman amplitudes. This is unlike the U(l) 
case, where the photon acquires a mass via an intrinsic Higgs mechanism. 
This has led to the speculations that QCDi may in fact exist in two phases 
associated with the weak and strong coupling regimes. In this picture, the 
large N limit would correspond to the weak-coupling limit ('t Hooft's phase), 
with massless gluons and a mesonic spectrum described by a Regge trajectory. 
In such a case, the Regge behaviour of the mesonic spectrum is compatible 
with confinement. In the strong coupling regime (Higgs phase), on the other 
hand, the gluons would be massive, and the original 5C/(A^) c -symmetry would 
be broken down to the maximal abelian subgroup (torus) of SU(N) C . 

The behaviour of the theory in the strong or weak coupling limits is rather 
subtle. The theory is asymptotically free. In the strong coupling limit, it is ex- 
pected to be in the confining phase: in the infinite infrared cut-off limit quarks 
disappear from the spectrum, which consists of mesons lying approximately 
on a Regge trajectory. 

The problem of screening and confinement can however be analysed along 
the lines of the U(l) case. However, unlike the U{\) case, the screening phase 
prevails in the non-abelian theory j|, 9. 
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1.1 QCD2 in the local decoupled formulation and BRST 
constraints 

The partition function of two-dimensional QCD in the fermionic formulation 
(before gauge fixing) is given by the expression 



J VA+VA^ J VipVtj;e iS[A ^ ] (1) 



with the action 



S\A,ip,i] = J d 2 x 



-trF^F^ + 4(id + + eA+^i + 4(id„ + eA-)i> 2 

(2) 

We can obtain a bosonic formulation of the theory, in such a way that struc- 
tural relations, hidden in the fermionic formulation are made clearer in the 
bosonic counterpart. We first make some useful change of variables, obtaining 
a formulation in terms of matrix-valued fields which decouple at the partition 
function level, but which are not totally decoupled, due to the gauge symme- 
tries of the theory. 

We can also arrive at a factorized form of the partition function ([l]) by 
parametrizing A± as 

eA + = U~ x id+U , A_ = Vid-V- 1 (3) 

as well as performing a chiral rotation, 

^ _, ^(0) = Urj3Xt ^ _, ^(°) = y-ifo . ( 4 ) 

We arrive at 

Z = ZP J VUVWJ G [W]J F [W]e lSYM[w] , (5) 

where Zp^ is the partition function of free fermions, and Sym is the Yang- 
Mills action given by 

S YM [W] = J d 2 x[E 2 + Ed+iWid-W' 1 )] (6) 

with W = UV. 

The field strength tensor Fo\ is given in terms of W and U or V in either 
of the two alternative forms 

F i = -^U- l [d + {Wid^W~ l )]U = -Vld-iW'hd+W^V- 1 . (7) 
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The Jacobian Jp is given, following Fujikawa, by 

J F [UV] = e- ir[uv] (8) 

while the determinant of the adjoint Dirac operator is 

J G [UV] = e- lCvr ^(detid + ) adj (detid.) adj , (9) 

where cy is the second Casimir of the group in question with the normalization 
facdfbcd = cySab of the structure constants, and r[C/V] is the Wess-Zumino- 
Witten action. Representing (det id±) a dj in terms of ghosts and choosing the 
gauge U = 1, we obtain 

z = zfz^z w , (10) 



where Z$ = and 



--J VWe~ i{1+Cv)r[w] e iSYM[w] = J VWe iSe ff [w] , (11) 

(12) 

where we introduced the W-effective action 

S eff [W]=S Y M-(cv + l)r[W} . (13) 

The effective action has two terms, one corresponding to a WZW action with 
a negative coefficient, and the Yang-Mills action written in the form (ph. 

We refer to (|l0|) as the "local decoupled" partition function. As seen from 
(JTT|) , the ghosts b±* are canonically conjugate to c±/ and have Grassmann 



parity +1. We assign to them the ghost number ghjfc = —1 and ghjfc = +1, 
respectively. 

The dimensionality of the direct-product space Tip <8>7^j <g) Hw associated 
with the partition function (|lC|) is larger than that of the physical Hilbert 
space of the original fermionic formulation. Hence there must exist constrains 
imposing restrictions on the representations which are allowed in 7~i p hy S - I n 
order to discover these constraints we observe that the partition function is 
separately invariant under the following nilpotent transformations [10, [H]] 



(14) 



W5W~ 


1 = -cL 0) , 


w^sw = 


- -c (0) 




5^ = 


cL°M 0) , = o , 


«^ 0) = 


x /(°) 

, 0W2 = 


- cV 0) 4 0) , 


6c^ = 


i{ c L 0) ,cL 0) },^ = o, 


fc<? = , 


5cf = 


1 fJO) JO) 

2 I + ' + 


5b W = 




5b {0) = , 


5bf = 
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where 5 denotes the variation graded with respect to Grassmann parity, and 
flzp are given by 



n_=-—x>_(w)(d + (Wid-W- 1 )) - (i + c v ) j-(W)+j. 
n + =-^v + (w)(d-(w-Hd + w)) - (i + c v ) J+(w)+j, 



with 



J_(W) = -l-Wid-W- 1 , J+(W) = —W'hd+W 

47T 47T 



3- 



(15) 



(16) 
(17) 



These transformations are easily derived by departing from the Yang-Mills 
action 0, [0J. 

The corresponding BRST currents, as obtained via the usual Noether con- 
struction, are found to be 



J T = tr c 



(0) 



n 



i 



(18) 



with d + J- = 0, and cL J + = 0. 

Remarkably enough, the nilpotent symmetries lead to currents J_ and J+ 
which only depend on the variable x~ and x + , respectively. 

The on-shell nilpotency of the corresponding conserved charges 

Q± = J dx l J±{x ± ) (19) 
follows from the first-class character of the operators fi± = tr (t a £l±). 



1.2 QCD2 in the non-local decoupled formulation and BRST 
constraints 



The partition function represented by the standard expression (|11| ) contains 
fields which are mixtures of massive and massless modes, of positive and nega- 
tive norms respectively, coupled by the constraints. In the following we disso- 
ciate these degrees of freedom by a suitable transformation. We shall thereby 
be lead to an alternative nonlocal representation of the partition function, 
useful for learning certain structural properties. 

Following Ref. [12], we make in (1(1) the change of variables E — > j3 defined 



by 

2tt 



d + E = ( i±£^ ) p-Hd+p . (20) 
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The Jacobian associated with this change of variables is 

VE = detiV + ((3)V(3 



(21) 



where we have suppressed the constant det d+ which will not play any role 
in the discussion to follow. Making use of the determinant of the fermionic 
operator in the fundamental or in the adjoint representation and representing 
(detz<9+) as a functional integral over ghost fields 6_ and c_, we have, after 
decoupling the ghosts, 

Z=zfz$£fl JVW J Vf3exp{-i(l + c v )[T[W] + T\0\ 

triP^d+pWd-W- 1 )}} (22) 



1 

47T 



x exp(ir[/3]) exp < i 



l + c v 
2vr 



1 

— tr 
2 



where 



Mo) 

^gh- 



(23) 



Using the Polyakov-Wiegmann identity, J|] and making the change of variable 
W -> PW = W, we are left with 



7 _ jr(0) jr(°) 7-. 7 

F ^gh ^gh-^W^P ' 



where 
and 



1 + cy 
2vr 



-tr 



VWexp[-i{l + c v )T{W]] 
is the partition function of a WZW field of level —(1 + cy). 



(24) 

(25) 
(26) 



1.3 Massive two-dimensional QCD 

The BRST symmetries of the physical states in massless QCD2 are also the 
symmetries which should be imposed on the physical states in the massive 
case. For massive fermions the functional determinant of the Dirac operator, 
an essential ingredient for arriving at the bosonised form of the QCD2 partition 
function, can no longer be computed in closed form, and one must resort to 
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the so-called adiabatic principle of form invariance[Q. Equivalently, one can 
start with a perturbative expansion in powers of the mass, as given by 



d xipip 



(27) 



Afterwards, we use the (massless) bosonization formulae and re-exponentiate 
the result. In this approach, the mass term is given in terms of a bosonic field 
g^p of the massless theory by[[b| 



S m = -MJM> = M V J tr(^ + g^) , 

where \i is an arbitrary massive parameter whose value depends on the renor- 
malization prescription for the mass operator. || 

Defining m 2 = M[i, we re-exponentiate the mass term. Going through the 
changes of variables leading to one arrives at the following alternative 

forms for the mass term when expressed in terms of the fields of the non-local 
formulation 

S m = m 2 j ti{gt- l p + f3~ l tg~ l ) . (28) 

The corresponding effective action of massive QCD2 in the non-local for- 
mulation reads 

S = S Y M[(3,B} + S m [g,f3,t} + T[g}+T[P}-(c v + l)T[t} + S gh + S gh . . (29) 

We thus see that the associated partition function no longer factorizes. 
Nevertheless, there still exist BRST currents which are either right- or left- 
moving, just as in the massless case. 

The action (|2^) exhibits various symmetries of the BRST type; however, 
not all of them lead to nilpotent charges. The variations are graded with 
respect to Grassmann number. The equations of motion obtained from action 
flU) read 

l-d+igd-g- 1 ) =m 2 (gfl~ 1 f3 - P^Eg' 1 ) , (30) 

-°- }L ^r d +^ d - t ' 1 ) =m 2 (tg- l r 1 ~ faZ- 1 ) , (31) 
-U+OScL/T 1 ) + iXd+iPB/3- 1 ) =m 2 (Pg£- 1 - S^ 1 /?" 1 ) , (32) 
-^d^(p- 1 d+f3)+iX\p- 1 d + l3,B] + 

47T 

i\d + B =m 2 (gt- l l3 - p^tg- 1 ) , (33) 

d\B =X(p- 1 id + P) , (34) 

d ± b T =0 , d±c T =0 , (35) 
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with an analogous set of equations involving a so called prime sector [[H]], and 
where A = Cy n +1 e. Notice that the mass term can be transformed from one 

Z7T 



equation to another, by a suitable conjugation. Making use of Eqs. (30-35), 
the Noether currents are constructed in the standard fashion: we make a 
general BRST variation of the action, without using the equations of motion, 
and equate the result to the on-shell variation, taking into account terms 
arising from partial integrations. The only subtlety in this procedure concerns 
the WZW term, which only contributes off-shell to the variation. The four 
conserved Noether currents are found to be 

J± =tr ^c±0+ - ^b±{c±,c±}^j , (36) 
J± =tr (c±n± - ^b±{c±,c±}^ , (37) 
where the constraints above generally denoted as O are given by 

Q_=(— gid^g- 1 - ^Ifid-t'- 1 + {6_,c_}^) , (39) 

V 47T 47T / 

&- = ( -l-Pid^p- 1 - CV A + 1 tid-t- 1 - \f3Bl3- 1 + {b-.c,}) , (40) 
\ An 4tt J 

n + =(^P'-Ud + P'-^^t'- 1 i d + t'-\p'- 1 B'P' + {b + ,c + }^ . (41) 



(see |10| , 11 1 for definitions concerning the primed sector, connected to the 
unprimed one by a nonlocal transformation). 
From the current conservation laws 

d T J± = , d±J T =0 , (42) 

one infers that f2_ and f2_ are right-moving, while + and + are left-moving. 
Indeed, making use of the equations of motion (|30|-35) one readily checks that 
the operators Q±, £l± satisfy 

d T fl ± = , d±&T = , (43) 

consistent with the conservation laws (|42|) . In Ref. |l(| it has been argued 
that gauge invariant bilinears are the physical states of the theory. There is 
of course the problem of whether we have to consider or not that the physical 
states are annihilated by the non local constraints, but that problem goes 
beyond the scope of the present paper. 



8 



1.4 Screening in two-dimensional QCD 

Let us reconsider the problem of screening and confinement. We shall concen- 
trate on the case of single flavour QCD, and merely comment on the general 
case at the end of the section. 

We proceed by first considering the case of massless fermions and compute 
the inter-quark potential. We introduce a pair of classical colour charges of 
strength q = q a t a separated by a distance L. Such a pair is introduced in the 
action (p9|) by means of the substitution 



i(JT l d + f3) a — » i{f3- l d + f3) a - (s(x - |) - 5{x + , (44) 

where a is a definite colour index. This adds the following term to the action^] 

V{L) = AS = S g -S = -(c v + l)q a (V(L/2) - B a (-L/2)^j . (45) 

The equation of motion for B a is now replaced by 

d\B a = i\{(i- l d + f5) a - (c v + l)q a (s(x - |) - 5{x + |)Y (46) 

which implies, upon substitution into the equation of motion for the /3-field, 

d + (^-^d-d+B + [d+B, B] + iAS 

^a_ + (c V + l)[ff > B])[J(x-^) + • (47) 

We look for solutions of ( |47D with a fixed global orientation in colour space.^] 
We thus make B a = q a f(x). This renders the problem abelian. We thus infer 
that the potential ( f45| ) has the form 

ViL) = (C " + 1) ^ ^(1 - e - 2 ^ } (48) 

which implies that the system is in a screening phase. 

We now turn to the case of massive fermions. Taking the external charge 
to lie in the direction t 2 of SU (N) space, our Ansatz for B a leads one to look 
for solutions with g, (3 and E parametrized as 

5 = e i2V?^r 2) p = e i2V?E<j 2j E = e -*2V5^a _ (49) 



1 This corresponds to minus the same term added to the Hamiltonian. 

2 Note that this is a non-trivial input, since we have no longer the freedom of choosing a gauge 
in which such an Ansatz could be realized. 
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The equations of motion (|30|-34) are replaced by a set of coupled sine-Gordon 
type equations. As is well known [^], the solution of the classical equations of 
motion in the bosonic version contains quantum mechanical information from 
the fermionic theory. We find, after solving them, the result || 



(cy + l?q 




(50) 



where 



cy 2yjir(cv + 1)<7 



m 



(cy + 1) ' ^ (1 + ea 2 ) ' 

8-Km 2 2- v /7rea(cy + l)q 

m+ — 16em 2 ' (1 + ea 2 ) 



2tt[(a 2 + (1 + e)2m 2 ) ± y(A 2 + (1 + e)2m 2 ) 2 - 8eA 2 m 2 ], 



Q±=q± 



G(x-^)-G(x+^ 



(51) 



(52) 



Thus we find two mass scales given by m+ and m_. Both these scales 
correspond to screening-type contributions if cy ^ 0. 

Next, we compare the above results for the potential with those obtained 
for the Schwinger model. In the abelian case, the combination of the matter 
boson ip and the negative metric scalar rj gives rise to the #-angle. That is, 
the combination $ = ip + rj = 8 appears in the mass term. When fermions 
are massless, the electric field and the matter boson decouple. However, due 
to a Higgs mechanism, the electric field acquires a mass and, therefore, a 
long-range force does not exist. This leads to a pure screening potential. On 
the other hand, for massive fermions, the electric field couples to the matter 
boson Yet, $ = ip Cv= o, and hence, it remains massless. The coupling 
to $ via the mass term is the origin of the long-range force (linearly rising 
potential) in the massive U(l) case, where the potential is confining. On the 
other hand, the expression ( |50|) for the potential indicates the absence of a 
long-range force in the non-abelian case. 

The abelian potential can also be obtained from ( |50|) by taking the limit 
cy — > 0. In this limit, the mass scale m_ tends to zero and we recover the 
linearly rising potential, signaling confinement. 
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1.5 Equations of motion and higher conservation equations 



We are now going to deal with the action S e f f given in ([11]) , obtaining further 
important information. 

Due to the presence of higher derivatives in that action, it is convenient 
to introduce an auxiliary field and rewrite it in the equivalent form (||). 

The equation of motion of the W^-field is easily computed. The WZW 



contribution has been obtained in [14], and the Yang-Mills action leads to an 
extra term. We obtain 

^r°* + ^1 {w-w-') - j^lwe-w-\ S+ (wa-w-*)] s 

We can list the relevant field operators appearing in the definition of the 
conservation law (|53|), that is 



,5^ + i),5 + ^ + ^5-^5,<M5] , (54) 

with d + lYf = 0, and J_ = WdJW~ x . It is straightforward to compute the 
Poisson algebra, using the canonical formalism, which in the bosonic formu- 
lation includes quantum corrections. We have 

{lj[(t,x), Ifat^l^u^-y 1 ) - a^S'ix 1 - y 1 ) , 
{l$(t,x), J^ / (t,y)}=J^.^ ] 5(x 1 - y 1 ) + 25 a 8 kj 5'(x 1 - y l ) , 
{jW^xlJ^y 1 )}^ , (55) 

where a = ^{cv + 1)> and the indices in the right hand side have been 
appropriately antisimetrized, as denoted by the bracketts in the indices of the 
current I and of the Kronecker delta <5y. We thus obtain a current algebra 
for acting on with a central extension. 

1.6 Dual case-non local formulation 

At the Lagrangian level, we find the Euler-Lagrange equations for f3 from the 
perturbed WZW action (^5[), that is, 



p- L 5(3 , (56) 



5A(P)=2(d^ 1 (p- 1 d + p)-[df(P~ 1 d + P),(P~ 1 d + P)])p- 1 5p . (57) 
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We define the current components 



J^p-^+P , (58) 
/=-4V^ 2 jf = -V^V" 1 ^/^) , (59) 

which summarize the /3 equation of motion as a zero-curvature condition given 

by 

[£, C] = [d + + 4, 5_ + f\ = 8^4- d+f + J0]=O . (60) 

This is not a Lax pair, as e.g. in the usual non- linear a- models, where J@ is a 
conserved current and a conserved non-local charge is obtained. However, to 
a certain extent, the situation is simpler in the present case, due to the rather 
unusual form of the currents. This permits us to write the commutator as a 
total derivative, in such a way that in terms of the current J_ we have 

d+ (4^ 2 / + d+d^ jt + [jf , d + f\) = o . (6i) 

Therefore the quantity 

f ( x - ) = 4TTfi 2 .f(x + ,x-) + Jt (x + ,x-) + [jP (x+ ,x-),d+J?.(x + ,x-)] 

(62) 

does not depend on x + , and it is a simple matter to derive an infinite number 
of conservation laws from the above. 

Canonical quantization proceeds straightforwardly, and as a consequence 
we can compute the algebra of conserved currents, which is analogous to (p5|). 

We are thus led to speculate whether two-dimensional QCD contains an 
integrable system^. The theory corresponds to an off-critical perturbation of 
the WZW action. If we write (3 = e*^ ~ 1 + i<p, we verify that the perturbing 
term corresponds to a mass term for <f>. The next natural step is to obtain the 
algebra obeyed by (|62[), and its representation. However, there is a difficulty 
presented by the non-locality of the perturbation. 

1.7 Algebraic aspects of QCD2 and integrability 

We saw that two-dimensional QCD, although not exactly soluble, in terms 
of free fields, is a theory from which some valuable results may be obtained. 
The 1/N expansion reveals a simple spectrum valid for weak coupling, while 
the strong coupling offers the possibility of understanding the baryon as a 
generalized sine-Gordon soliton. 
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All such results point to a relatively simple structure, which could be 
mirrored by an underlying spectrum generating algebra. In fact such an alge- 
bra does exist. In the above-mentioned case of the large- N expansion of pure 
QCD2, one finds a spectrum generating algebra related to area-preserving 
diffeomorphisms of the Nambu-Goto action. For gauge-invariant bilinears in 
the Fermi fields this algebra can be constructed [15|. It also appears in the 
description of the quantum Hall effect. Moreover, pure QCD2 is equivalent 
to the c = 1 matrix model, which also has a representation in terms of non- 
relativistic fermions. The problem is also related to the Calogero-Sutherland 
models. The mass eigenstates build a representation of the algebra as 
found in IB]. 

In spite of the hints toward a possible integrable structure in two-dimen- 
sional QCD the problem still remains largely openflq]. A possible explanation 



was also given in [17|. The problem points to a better understanding of the 
mesonic spectrum of the theory, 't Hooft's results implying an infinite number 
of mesons obeying Regge behaviour is compatible with confinement. However, 
results obtained by several authors [||, ^, |l^] point to a screening phase. In that 
case the quark anti-quark potential looks too shallow to support an infinite 
number of bound states. 

Here we investigate the decay amplitudes of two-dimensional QCD using 
a refined numerical analysis, which will permit to account for the details of 
the dynamics. 



2 Numerical backup 

We have investigated the behaviour of the decay amplitudes of 't Hooft mesons 
in the theory, in the large N limit. 

One of the strong reasons to study decay amplitudes is to learn about the 
stability of the 't Hooft mesons, namely answer whether they might repre- 
sent soliton solutions of two-dimensional QCD. Moreover, these mesons are a 
probe into the long range force: indeed, if two-dimensional QCD is confining 
there is a linearly rising potential which acomodates an infinite number of 
bound states. On the other hand, screening would imply a shallow potential, 
flatening at infinity, implying a finite number of mesons. 

Decay amplitudes in the framework of perturbation theory in the inverse 
number of colours has been studied by a few authors |l9[ |. The essential 
ingredient is first the solution of 't Hooft equation in order to find the bound 
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Fig. 1: Amplitudes for the decay processes n — > 1 + 1 {n = 5, 8, 9, . . . , 18) 
versus the squared fermion mass m 2 , varying from to 0.1. 



state wave functions and meson mass-square eigenvalues /x 2 : 

^4>{x) = -^-Mx) - P f 1 dy-^-, , (63) 



x(l-x) Jo (x-y) 

where P denotes de Cauchy principal value prescription. Here 7 refers to 
a square mass scale corresponding to the quark-antiquark pair with equal 
masses m 2 , where the massless limit is obtained for 7 = @. 

This integro differential equation cannot be solved by analytic methods. 
However, using a sensible numerical method presented in |2(| one can arrive at 
quite reasonable results for the eigenvalues and wave functions. This method 
makes use of a spline procedure to work out the exact wave functions on an 
adaptive grid. The numerical representation of the eigenfunctions are used 
to evaluate the decay amplitudes from an initial meson state n into two final 
states ni and ri2- In the leading and next to leading order in the 1/N expansion 



they are given by [EU 22 



I [ w X x — VJ 

(n|ni,n 2 )=(l - C)- / dx<j> n (x)(j> ni (—)$ n2 {- ) 

1 — w Jo w 1 — w 
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Fig. 2: Amplitudes for the decay processes n — > 1 + 1 (n = 5, 8, 9, . . . , 18) 
versus the squared fermion mass m 2 , varying from to 2.35. 



I f X — TV X 

-(l-C)-/ dx<f) n (x)(f> n2 (- )$m(-) 

W Jw 1 — W W 



1 i -i ^-,N| /ri2 /"'" j t t „\ j. I ' 



+-(1-C)-^/ cte0 n (x)0 m (-), (64) 
iv 1 — u> Jn w; 



o 

where C denotes de interchange of final states, 

" x ^dy) 
/o 3 (x - y) 



*»(*) = I dy^L- (65 ) 



and w is a kinematic parameter, whose on-shell values w + and u>_ correspond 
to the right moving and left moving of the final state n\. 



As observed in [17] the higher-order corrections for the amplitude are al- 



ways multiplied by the factor 

fn 2 = [ dx(f> n2 (x), (66) 

where <f) n2 (x) is the 't Hooft's eigenfunction corresponding to the decaying 
state 7i2, which vanishes for massless fermions (for n,2 ^ 0). 
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Fig. 3: As in figure 1, for the decay processes n — > 1 + 2 (n = 8, 9, . . . , 18). 



With the knowledge of the eigenfunctions the numerical calculation of 
these integrals is straightforward and we obtain the amplitudes in function 
of the outgoing momenta w. From here on we report the behaviour of these 
amplitudes for the on-shell parameter to + , 

A(w + ) = {n\n 1 ,n 2 ) w = w+ . (67) 

For this end we had to use an interpolation and extrapolation numerical al- 
gorithm |23|] to obtain a precise value of K(w+) from amplitudes evaluated at 
a discrete set of values for w. 

First we compute the decay amplitude of a level- n state into two level- 1 
states. Notice that level-0 state decouples for zero fermion mass (s ee Jl7|) and 
at large n we also expect these decay amplitudes to approach zero |17|> |j] . As 
argued in [17] we expect a small anomalous amplitude for vanishing mass, and 



a nonvanishing result for the massive fermion case. Moreover, there might be 
a new physical description at some nonvanishing value of the fermion mass as 
compared to the (dimensionfull) coupling constant, namely m 2 ~ — , since 't 
Hooft's solution present a quark propagator for m 2 < — . In Fig. 1 we show 
the results obtained for quark-antiquark pairs with masses m 2 ranging from 
to 0.1, while in Fig. 2 we extend this range up to 2.35. Interestingly enough 
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Fig. 4: As in figure 2, for the decay processes n — > 1 + 2 (n = 8, 9, . . . , 18). 



are the several zeroes of the amplitudes for a fermion mass comparable to 
the coupling constant e, but they do not coincide for different decays. This 
behaviour motivated us to continue exploring the decay amplitudes involving 
high excited levels n (n = 8,9, ...,18) into the next low-lying levels n\ = 
l,7i2 = 2 (Fig. 3 and 4) and, n\ = 2,ri2 = 2 (Fig. 5 and 6). We observe 
that the feature described above is reproduced in figures 4 e 6. Whether these 
zeroes signalize any new physics hidden by anomalies is an issue about which 
we can unfortunately only speculate. Furthermore, as the bosonic state level 
increases (resp. decreases), it becomes stable at higher (resp. lower) mass 
values due to the lack of possibility in terms of kinematical variables, to fulfill 
the energy momentum conditions. This can be seen as interrupted lines in 
figures 2, 4 and 6. It is a question under investigation whether we obtain 
further stability conditions while correcting the mesonic wavefunctions and 
spectrum in order to cope with the screening mechanism, and not using a 
confinement potential between quarks and antiquarks, as implicit in 't Hooft's 
solution. 

It is worth to mention the decays n — > 2 + 2, namely an arbitrary state 
n into two mesons at the second excited state present these issues in a more 
enhanced way. Among the studied states, the first to decay according to this 
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Fig. 5: As in fig. I, for the decay processes n — > 2 + 2 (n = 8, 9, . . . , 18). 



mode is the eighth excited state. The smallness of the decay amphtudes for 
m ~ can aiso be verified and, for not very big vaiues of n, it is clear iy seen 
how the sequence of forbidden decays due to the violation of the kinematic 
condition appears as the mass increases. 

3 Conclusions 

After more than two decades of study, two-dimensional QCD largely remains 
a challenge for the complete understading of gauge theories. Nevertheless very 
deep concepts emerged and a plethora of information concerning its dynamical 
understanding could be gathered. 

't Hooft's solution taught us how to perform the 1/N expansion of the 
theory, but also its limitations. From that solution emerged a massless gauge 
field, contrary to Schwinger model inputations, and a spectrum compatible 
with confinement of the fundamental degrees of freedom, namely the quarks. 
The 1/N expansion has been further studied, and here in particular we derived 
the numerical results of the decay amplitudes. 

It has been revealed that the structure of these amplitudes is very simple, 
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Fig. 6: As in figure 2, for the decay processes n — > 2 + 2 (n = 8, 9, . . . , 18). 



and in most of the cases rather small, especially for the massless case, thus 



compatible with the speculations set forward in Ref. [17]. Indeed, we see that 
in figures 1 through 4 all amplitudes are small for zero mass, but non vanishing, 
raising the suspicion of a higher conservation law, but with the presence of 
a small quantum anomaly, effective for the decay of low levels of the meson 
mass. Therefore, it turns out to be very important to check whether the 't 
Hooft hierarchy is compatible with the screening rather than the confining 
potential, in which case the infinite number of states are substituted by a 
finite number of mesons. It is to be checked in such a case whether decays 
might be further supressed. 

Thus, the detailed confirmation of the Regge behaviour and higher states 
in the framework of the 1/N expansion, which is seemingly incompatible with 
the screening potential obtained by several independent authors || || [l^, 18, 



25, 26] requires a better understanding of the spectrum, and leads to the 



suspicion that it contains a finite number of mesonic states. Moreover, it has 
been recently pointed out that solitons do not exist in the model, strenghening 
the conclusions drawn in [17] that the symmetries suspected in [12] where 
anomalously realized. Following this vein, we confirm that it is essencial 
to obtain an alternative means of checking the spectrum of the theory and 
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comparing to 't Hooft's result. In case we arrive at a confirmation of the 
latter, we need to explain the dual description of screening, as obtained from 
computation of Wilson loops of semiclassical potential and confinement from 
the Regge behaviour of the mesonic states. 

From the results of [jlO| we also conclude that physical states of the the- 
ory can be constructed out of bilinears. This points to a confirmation of 't 
Hooft's spectrum. On the other hand, as we have seen in the introduction this 
depends on the rather unclear structure of the non local constraints we have 
obtained, which do not inspire confidence due to the lack of mathematically 
sound theorems concerning them. 

Further analysis of the diagrams such as figures 5 and 6 point also to a 
simplification of the amplitudes for large values of the mass, and apparently 
no transition or anomalous behaviour for the fermion mass comparable to the 
coupling constant. At that point, a transition between weak and strong cou- 
pling has been long suspected^], as signalized by the existence of a tachyon 
pole in the (formal) quark propagator. This points once more to an urgent 
necessity of further analysis of the large N limit of two dimensional QCD. A 
point missed in the large iV analysis is the fact that for finite iV the gauge bo- 
son has a mass generated by the Higgs mechanism, analogous to the Schwinger 
model case, and easily derived from the pseudo divergent of the Maxwell equa- 
tion with use of the anomaly equation [0], but the ensuing mass vanishes for 
large N. 

Although conclusions are still insufficient to have a complete physical de- 
scription of the physics of the model, especially in view of the efforts spent 
in the problem, progress have been achieved, and in fact the results already 
spilt off to the three-dimensional case [27], where the screening mechanism 
agains prevails, leading to the suspicion that it is physically more appealing 
to the theory to form kinkstates to take care of the long range force as pro- 
posed by p8| rather than to leave for the naked gauge field to pull the quarks 
together by means of a long range force. If this kind of mechanism works also 
in the four dimensional case the physical consequences will be far reaching, 
especially concerning the baryonic and mesonic spectrum. 
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